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ABSTRACT 

Using radio data at 1.4 GHz from the ATCA we identify five head-tail (HT) 
galaxies in the central region of the Horologium-Reticulum Supercluster 
(HRS). Physical parameters of the HT galaxies were determined along with 
substructure in the HRS to probe the relationship between environment 
and radio properties. Using a density enhancement technique applied to 582 
spectroscopic measurements in the 2° x 2° region about A3125/A3128, we 
find all five HT galaxies reside in regions of extremely high density (>100 
galaxies/Mpc"^). In fact, the environments surrounding HT galaxies are 
statistically denser than those environments surrounding non-HT galaxies 
and among the densest environments in a cluster. Additionally, the HT 
galaxies are found in regions of enhanced X-ray emission and we show that 
the enhanced density continues out to substructure groups of 10 members. 
We propose that it is the high densities that allow ram pressure to bend 
the HT galaxies as opposed to previously proposed mechanisms relying on 
exceptionally high peculiar velocities. 

Key words: radio continuum: galaxies; galaxies: clusters: general; galax- 
ies: clusters: individual: A3125, A3128 



1 INTRODUCTION 

The term "radio galaxy" is usually used to describe a 
class of objects which have characteristic jets of highly 
relativistic particles blasted out from the host galaxy. 
Although most if not all galaxies have some radio 
emission, the term "radio galaxy" is usually reserved 
for those galaxies that h ave a high radio-to-optical lu- 
minosity ratio ( Kellermann et al.|[l989. '). 

Head-tail (HT) galaxies are a sub-class of radio 
galaxy whose radio jets are distorted such that they 
appear to bend in a common direction. This gives rise 
to the so-called "head-tail" morphology where the ra- 
dio jets are bent back to resemble a tail with the lu- 
minous galaxy as the "head". (Th e term "head-tail' 
was coined by iMilev et all ll 19721) but the first H T 
galaxy was discovered by Rvle fc WindramI (llQGSh .) 
Although considered somewhat exotic in their early 
years, HT galaxies are now considered quite common. 

HT galaxies are generally detected in dynami- 
cal, non-relaxed clusters and are preferentially found 
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towards the centre of the cluster potential well 
l|Klamer et all |2004 ) or in regions of enhanced X- 
ray emission indicativ e of a deep potential well 
(|jones fc FormanI Il984h . Although it is common to 
find HT galaxies in clusters, the number of HT galax- 
ies per cluster is low with most unrelaxed clusters that 
have been probed yielding only one or two HT galaxies 
per cluster. 

HT galaxies are g enerally FRI in morphol- 
ogy l|Venturi et al.l 1 19981 ) which, assuming that the 
FRI/FRII break is caused by environment, is consis- 
tent with HT galaxies being found in clusters as their 
denser e nvironments effectively restrict the outfiow of 
the jets (|Owen fc Ledlowlll994h . 

The bent morphology of jets in HT galaxies has 
always been attributed to one of two possible mecha- 
nisms relating the relative velocity of the host galaxy 
to the intra-cluster medium (ICM). In the first mech- 
anism the host galaxy has a larger than expected pe- 
culiar velocity, leading to ram pressure on the jets 
l|Milev et al.ll 19721 : iRudnick fc Owenll976l ). where ram 
pressure. Pram, is defined as: 
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Pram OC pV , (1) 

where p is the density of the ICM, and v is the 
relative velocity be tween the galaxy and the ICM 
(|Gunn fc Gotdl 19721 ). 

In the second mechanism, strong winds in the 
ICM, believed to be caused by dynamical interac- 
tions such as cluster-cluster mergers or accretion 
of groups on to clusters, bend the jets. The phe- 
nomenon is often referred to as "cluster weather" 
(|Burnsl[l99i ). These mechanisms have been invoked 
to explain the difference in the degree to which 
the jets are bent, with Narrow- Angle Tailed galax- 
ies (NATs) believed to be formed via ram pressure 
in cases where the host galaxy has a peculiar veloc- 
ity of the order of 600 km s~^ and is em bedded in a 
hot, dense plasma (Vcnkatosa n et a"l]|l994l '). and Wide- 
Angle Tailed galaxies ( WATs) believed to be the result 
of "cluster weather" jKlamer et al.ll2004 ). NATs are 
most likely to be ca used by th e passage of the galaxy 
through the cluster (|Odea fc O wen 1985) while WATs 
are more likely to be associated w ith dominant cluster 
galaxies (jOwen fc Rudnic^ll976h . 

Although the distinction between ram pressure 
and ICM effects, such as "cluster weather" (|Burnsl 
1 19981 ). is often drawn in literature, we note that clus- 
ter weather still implies ram pressure as the cause of 
the jet bending. The true distinction between the ram 
pressure theory and the cluster weather theory lies in 
the velocities of the HT galaxies relative to the kine- 
matic centre. Henceforth, when we refer to the ram 
pressure theory we are indicating that the galaxy has 
a high peculiar velocity and is in fact travelling at a 
velocity which differs significantly from the statistical 
distribution of velocities for other galaxies in the clus- 
ter. Cluster weather on the other hand will mean that 
the galaxy does not have a significantly different ve- 
locity, but rather it is the motion of the ICM that is 
causing the bent morphology. 

In both cases it is the relative motion of the host 
galaxy with respect to the ICM that causes the bend- 
ing. Consequently when only viewing the radio mor- 
phology it is not possible to distinguish between a 
galaxy with high peculiar velocity moving through an 
undynamical, stationary ICM, and a galaxy with the 
ICM rushing past. However, spectroscopic redshifts 
can be used to determine the distribution of veloci- 
ties for cluster galaxies and, via a statistical analysis, 
the relative importance of the host galaxy's motion 
with respect to its neighbours can be established. 

This paper explores the environment of a sam- 
ple of HT galaxies and examines their usefulness as 
a tool to investigate dynamical conditions in the host 
clusters. We examine the well-studied major merger 
A3125/A3128 which resides in the central region of 
the Horologium- Reticulum Supercluster (HRS). Us- 
ing a large spectroscopic dataset in conjunction with 
multi-frequency radio imaging we present a quanti- 
tative analysis of the environment of HT galaxies as 
compared to samples of both other cluster galaxies 
and other cluster radio galaxies. The paper is laid out 
as follows; Section 2 discusses the spectroscopic and 



radio data in the region surrounding A3125/A3128 
and the parameters associated with the detected HT 
galaxies; Section 3 presents the substructure detec- 
tion; Section 4 gives the results of our statistical anal- 
ysis of the environments surrounding HT galaxies and 
Sections 5 and 6 present the discussion and conclu- 
sions. 

This paper uses Ho = 70 km s~^ Mpc~^ and JIm 
= 0.27. This gives a scale of 1.191 kpc/arcsecond at 
18500 km s"\ the mean redshift of the HRS. 



2 HEAD-TAIL GALAXIES: DETECTION 
AND MORPHOLOGY 

To investigate the relationship between HT galax- 
ies and the dynamics of the cluster, we study HT 
galaxies in A3125/A3128, a major merger in the cen- 
tral region of the Horologium-Reticulum supercluster 
(HRS). The HRS is the second largest supercluster in 
the local universe (out to 300 Mpc), preceded only 
by the Shapley supercluster. The cores of rich super- 
clusters are ideal environments to study cluster merg- 
ers because the high peculiar velocities induced by 
the enhanced local dens ity favours cluster-c luster and 
cluster-group collisions (jVenturi et ahlfioOOl ). 



2.1 Spectroscopic Data 

We compiled a catalogue of 2181 spectroscopic red- 
shifts in the central region of the HRS. Data 
were obtained from all previou sly published surveys 
l|Rose et al.ll2002l : iFleenoillioO^ ) in addition to an un- 
published catalogue of spectroscopic data obtained at 
ANU by T. Mathams in 1989 and with data from 
the NASA Extragalactic Database (NED). A spatially 
and velocity limited subset of the 2181 galaxies cov- 
ering 2° X 2° in area and 15,000 km s ~^ < cz < 
22,000 km s~^ in velocity centred on the major merger 
A3125/A3128 yields 582 redshifts. 

When compiling a catalogue from several differ- 
ent sources which have overlapping observations of the 
same galaxies, care must be taken to ensure that in- 
dividual objects do not appear more than once in the 
final catalogue or that close, but independent objects 
are not incorrectly conflated. To avoid these problems 
we employ a visual inspection method for all cases 
where objects from different sources appear within the 
spatial and recessional uncertainties of the catalogues 
18" and 400 km s~^). In such cases postage stamps 
from the Digitized Sky Survey (DSS) are inspected. 
This is an important step as a completely automated 
compilation method might inadvertantly remove true 
close-proximity galaxies. 

The final catalogue of 582 galaxies has a magni- 
tude limit of mi,=23 and is ~90% complete to mf, = 18 
and ~60% complete to mi,=19.5. The optical coverage 
of the HRS region appears uniform (see Figure [l}. 

The resultant spectroscopic catalogue, when com- 
bined with the available radio images, provides an 
ideal dataset with which to investigate the environ- 
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Figure 1. Top left: Skyplot of the HRS region with cutoff = 18. Top right: Skyplot of the HRS region with m;, cutoff = 
18.5. Bottom left: Skyplot of the HRS region with mi, cutoff = 19. Bottom right: Skyplot of the HRS region with m;, cutoff 
= 19.5. Contours are set at 10% intervals. The crosses mark the positions of the HT galaxies. 
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Figure 2. Histogram of galaxy velocities from our position and velocity limited catalogue (of 582 galaxies). The shaded 
histogram plots the 46 radio galaxies. The velocities of our HT galaxies are shown with dotted lines. It can be seen that 
the HT galaxies are clearly in the main component of A3125/A3128, as are the radio galaxies. 



ment of HT galaxies. Figure [5] shows the distribution 
of velocities of our catalogue. 

2.2 Radio Observations 

Radio images at 1.4 GHz of the central 3.28 square de- 
grees surrounding the A3125/A3128 system (centred 



at 03:29:45, -53:00:00) were obtained with the Aus- 
tralia Telescope Compact Array (ATCA). The obser- 
vation was a 40-pointing, linear mosaic consisted of 
2 X 128 MHz observations centred at 1344 and 1432 
MHz. The data were obtained in three different ar- 
ray configurations over 39 hours in 2001 (see Table 
[T]). The primary calibrator was PKS B1934-638 and 
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the phase calibrator was PKS B0302-623. The syn- 
thesized beam size is 9.90 by 7.87 arcseconds with 
a position angle of 31.9°. The mean RMS noise for 
the image (Figure [Sjl is 0.1 mjy beam~^ and lowest 
contours are set at 0.3 mJy beam~^, three times the 
RMS noise. Ov er 500 radio galaxies are detec ted in the 
1.4 GHz data l| Johnston- Hollitt et all 120081 ') of which 
46 are identified with optical galaxies having redshifts 
found in our spectroscopic catalogue that placed them 
in A3125/A3128. 

Low frequency data for the region were also avail- 
able from the Sydney University Molonglo Sky Survey 
(SUMSS) (Mauch ct al. 2003), which used the Molon- 
glo Observatory Synthesis Telescope (MOST) to as- 
semble a comprehensive picture of the Southern radio 
sky {5 < -30°) at 843 MHz. Data obtained with the 
MOST have a 3 MHz bandwidth and the resultant 
images from SUMSS have a resolution of the order of 
1'. 



2.3 Identifying HT Galaxies 

Five HT galaxies are identified in the 1.4 GHz 
data, including two that w ere previously detected 
(jjohnston-Hollitt et al.ll200^ ). Typically only one or 
two HT galaxies tend to be found in non-relaxed clus- 
ters in the current literature. However, the available 
radio surveys covering clusters tend either to be quite 
shallow, or have low sensitivity to e xtended emission 
(e.g. the NVSS (jOondon et al.lll998l )) and hence only 
the brighter extended sources are detected. In a com- 
parably shallow survey to those commonly used we 
would only be a ble to see the two HT galaxie s previ- 
ously studied bv lJohnston-Hollitt et all (|2004l ) as they 
are t he brightest H T galaxies of our sample. Interest- 
mg Iv. iMillerl (|2005h surveyed 7 square degrees in the 
central region of the Shapley Supercluster to a similar 
depth to our data (80 /xJy beam~^ at a resolution of 
16"), and also found five HT galaxies. 

The optical and radio parameters of the HT 
galaxies are discussed below and Table [2] summarises 
the respective frequency data. 



2.3.1 Optical Counterparts 

The positional data for the HT galaxies is taken as the 
coordinates of their optical counterparts; these val- 
ues were obtained from NED. The HT galaxies all lie 
at comparable distances and have similar magnitudes. 
The velocities of the HT galaxies are marked by dot- 
ted lines in Figure [S] 



2.3.2 Flux Densities 

We used the Karma (|Goochl package Kvis to 

determine the fiux densities of the HT galaxies at both 
available frequencies. The error in the flux density is 



Because the flux density of our radio galaxies is 
so high relative to the RMS noise, the uncertainty in 
total flux density is dominated by the 10% residual 
calibration error typical for the ATCA. 



2.3.3 Power and Fanaroff- Riley Classification 

Using the radio power formula from iHoggI l|l999l ) we 
calculate the radio powers for the five HT galaxies. 
Their powers range between 8.37 x 10^^ W Hz"^ and 
4.46 X 10^" W Hz-i. We find all of our HT galaxies 
to be either FRI or at the transition between FRI and 
FRIT, which is consiste nt given the abso lute magni- 
tudes of these galaxies l|Owen et al.ll2000l). This is in 
agreement with previous literature (eg. Venturi et al.l 
11998 ). 



3 SUBSTRUCTURE DETECTION 

Substructure, which may be defined as the presence 
of two or more clumps of matter within a cluster, is a 
clear sign of incomplete relaxation and is hence indica- 
tive of the dynamica l nature of a cluster of galaxies 
l|Pinknev et al.ll 19961 ). 

Substructure can be detected in two-dimensions 
via analysis of X-ray surface brightness data, but to 
investigate the full three-dimensional nature of a clus- 
ter a large spectroscopic survey is required. 

By identifying substructure in the cluster we can 
investigate the relationship between the location of the 
HT galaxies and the u n derlyin g dynamical structure. 

Caldwell fc Rose! l|l997l ) showed that the 



AS = \/{Sb)2 + a2 



(2) 



where S is the fiux density, b is the residual calibration 
error, and a is the RMS noise. 



A3125/A3128 system is in a post-merger state and 
that A3125, the smaller of the two clusters, is more 
dispersed due to its recent passage through A3128. 

Substructure in the cen tral region of the H RS has 
been studied previously by iRose et ah I l|2002l ). Using 
their spectroscopi c data for 532 galaxies in this region, 
iRose et al.l l|2002l ) identified four groups on the high 
velocity side of A3125/A3128. These groups were se- 
lected on the basis of their clumping in either RA-cz 
space or Dec-cz space. We note that in this context 
"groups" refer to substructure groups and not galaxy 
groups. 

Furthermore, iRose et al.l l|2002l ) noted the pres- 
ence of a void 4000 km s"'^ on either side of the clus- 
ters. They interpret this as a consequence of the large 
gravitational potential of the cluster and hence the 
voids are indicative of galaxies being absorbed from 
the peripheries into the clusters. 

The substru cture detection method used by 
iRose et all (|2002l ) is limited in that groups can only be 
oriented along the RA, Dec or cz axes. This means any 
"infalling groups" that are falling "diagonally into the 
cluster" , for example, will not be as easily detectable. 

We believe that rather than id entifying su b struc - 
ture in the A3125/A3128 region, iRose et al.1 (|2002| ) 
have merely identified parts of the cluster. For this 
re ason we have chos en not to repeat the method used 
bv lRose etU] l|2002l ) when detecting substructure. In- 
stead we use density enhancements in our catalogue 
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Table 1. Details of the 1.4 GHz ATCA observations. Columns 1 and 2 give the telescope configuration and observing 
date respectively. Column 3 is the integration time in minutes, column 4 is bandwidth in MHz, column 5 gives the central 
frequencies for each IF and column 6 is the uv-range covered in kilowavelengths. Observations were carried out by R. W. 
Hunstead, M. Johnston-H oUitt, J. A. Rose, W . Christiansen, & M. Fleenor. Full details of the observations and reduction 
are available in (iJohnston-HoUitt et al.ll2008l ). 



Configuration 


Date 


Integration 


Bandwidth 


Frequencies 


UV-range 






(min) 


(MHz) 


(MHz) 


(kA) 


6A 


8th & 10th December, 2001 


880 


2 X 128 


1344 & 1432 


0.50 - 20.72 


6D 


27th November, 2001 


840 


2 X 128 


1344 & 1432 


1.57 - 27.72 


1.5D 


19th & 20th November, 2001 
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2 X 128 


1344 & 1432 


0.36 - 27.43 
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Figure 3. Central region of the 3.28 square degrees 1.4 GHz image of A3125/A3128 overlaid on optical data from the DSS. 
The RMS of the image is 0.1 mjy beam~^ and the contour is at three times the RMS noise. The five HT galaxies identified 
in the image are denoted with boxes. We note that Galaxy A appears to have a companion galaxy nestled between its lobes. 
This galaxy is in fact in the background and hence not considered in this study. The X-ray peaks of A3125 (lower right 
with three peaks) and A3128 (upper left with two peaks) are marked with crosses - the size indicates the strength of the 
X-ray emission. 
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Table 2. Optical and radio parameters of the HT galaxies. Column 1 is the identifier of the HT galaxy while column 2 
gives its NED ID. Columns 3 and 4 give the Right Ascension and Declination in J2000 coordinates while columns 5, 6, 7 
and 8 give the redshift, luminosity distance (Wright 2006), apparent magnitude and absolute magnitude respectively. All 
magnitudes are b^, except for galaxy E which is B^. Columns 9 and 10 give the flux density and power of the galaxies at 
1.4 GHz respectively and column 11 gives the Fanaroff-Riley morphology class. 



A 


2MASX 


J03275206- 


-5326099 


03 


27 


52 


-53 


26 


10 


17745 


265 





16.20 


-20.92 


245 


± 


25 


1.84 


X 


1024 


FRI/II 


B 


2MASX 


.103272476- 


-5325179 


03 


27 


25 


-53 


25 


18 


18827 


282 





15.85 


-21.40 


529 


± 


53 


4.46 


X 




FRI/II 


C 


2MASX 


.103301366- 


-5237304 


03 


30 


14 


-52 


37 


28 


19417 


291 


2 


16.24 


-21.08 


35 


± 


4 


3.12 


X 


1024 


FRI 


D 


2MASX 


J03301522- 


5234124 


03 


30 


15 


-52 


34 


10 


17008 


253 


6 


16.94 


-20.08 


24 


± 


2 


1.62 


X 




FRI 


E 


IC 1942 


NEDOl 




03 


27 


54 


-52 


40 


35 


17964 


268. 


6 


17.27 


-19.87 


11 


± 


1 


8.37 


X 


1024 


FRI 



to determine the location of substructure in this re- 
gion. Using a physical parameter such as density en- 
hancement provides a more accurate analysis of the 
substructure within A3128/A3125. 

3.1 3D Density Enhancements 

To identify the location of substructure in 
A3125/A3128 we introduce our own method of 
identifying substructure using density enhancements. 

We searched for density enhancements in our cat- 
alogue of galaxies using the HIPAS S group- f inder, a 
group finding program written by IStevend l|2005l fl. 
This group finder needs only sky position and veloc- 
ity to operate, and finds groups using the hierarchical 
algorithm (|Tullv,.1987 '). 

Specifically , the group-finder uses the group- 
ing equations of lGourgoulhon. Chamaraux. fc Fougud 
([l9.92) (hereafter GCF92) to determine the likely 
three-dimensional separation of two galaxies us- 
ing their on-sky positions and recessional velocities. 
The GCF92 equations compute separation in three 
regimes: 

(i) When the velocities of the two galaxies differ by 
less than some user-specified parameter Vi, the veloc- 
ity difference is considered to be entirely peculiar, and 
thus the separation is calculated using only plane-of- 
sky information. 

(ii) When the velocities of the two galaxies differ 
by more than some user-specified parameter Vr, then 
the galaxies are considered to be ungroupable; ie. they 
are separated by Hubble expansion and are not likely 
to be bound to each other. 

(iii) When the velocities of the two galaxies differ 
by an amount between these two extremes, then the 
separation is calculated using a combination of plane- 
of-sky and velocity information. 

At the start of the procedure, each galaxy is con- 
sidered as an individual "unit" . The group-finder cal- 
culates the separation between every pair of units us- 
ing the GCF92 equations, and from this the radius of 
the groups (the average galaxy-galaxy separation of 
the group members) that would result from each pair. 

^ http://www-ra.phys.utas.edu.au/ jstevens/research.html 



The galaxy number density for each potential group is 
then calculated simply as the total number of galaxies 
in the two units divided by the volume. This differs 
from the method of GCF92, who calculated mass den- 
sity, but this is the only significant difference. 

The pair of units with the highest number den- 
sity are then grouped into a single unit, which re- 
places the two original units. The separations between 
each of the remaining pairs of units in the catalogue 
are then recalculated and the process begins again. It 
should be noted that a unit with more than one galaxy 
does not hide any information: when the radius of a 
group is calculated, the group-finder always considers 
the individual galaxies. The group-finder stops when 
the number density of the densest pair is below some 
user-specified stopping density. For each group that it 
identifies, the group-finder calculates both the radius 
as calculated from the GCF92 equations and the pro- 
jected radius, using just the plane-of-sky information. 

Because the group-finder uses number density to 
identify groups, it can be thought of as a type of 
nearest-neighbour finder. Indeed, if Vi is set to a suffi- 
ciently large value, and Vr was set to the same value, 
the group-finder would use only plane-of-sky informa- 
tion and would identify only those groups which ap- 
pear as dense regions on the sky. However setting Vi 
lower will allow the group-finder to identify structure 
along the line of sight. 

A comparison could be made to the non - 
parametric algorithm DEDICA (|Pisanil [l993 . 1 19961 ) . 
which was designed to identify clusters and determine 
whether they have substructure. However, the work 
presented in this paper is less concerned with the iden- 
tification of the clusters themselves, and more con- 
cerned with the environment immediately surround- 
ing each galaxy in the cluster. As the hierarchical al- 
gorithm gives information about clustering strength 
at all scales, it is particularly suited for our purposes. 

3.2 2D Density Enhancements 

To minimize the uncertainties that accumulate when 
using velocity as a pseudo-spatial parameter, we re- 
peat the density enhancement study using projected 
radii instead of 3D density. That is, the substructure 
groups are still selected using exactly the same method 
as above, but we measure the projected radii of the 
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Table 3. Cut-ofT densities for HT galaxies. 



Vi = 1000 km s" 
HT galaxy Cut-ofT density 
(galaxies/Mpc'^ ) 



A 
B 
C 
D 
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60680 
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20072 

64919 
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Figure 4. Results of the HIPASS group finder at V; = 1000 
km s-i. The X-axis gives the density in galaxies /Mpc^ on a 
logarithmic scale while the y-axis is the number of galaxies 
that can be found in groups. 

groups to determine relatively whether HT galaxies 
reside in different environments to non-HT galaxies. 

For each galaxy in the catalogue, as more galaxies 
are grouped with it, we can track how the projected 
radii increases as well, in effect showing us how density 
changes with increasing group membership. 



4 RESULTS 

From Figure[2]we can deduce that all the HT galaxies 
lie in the main peak of the velocity distribution. Thus 
it appears that the HT morphology is not a conse- 
quence of high peculiar velocity. 

Figure |4] shows the fraction of the spectroscopic 
catalogue assigned to a group as a function of stopping 
density. It is clear that at very low densities (~10~^ 
galaxies/Mpc^), all the galaxies in the catalogue are 
placed in groups and at very high densities (~10^ 
galaxies/Mpc^) almost no galaxies are in groups as 
expected. 

The cut-off densities of the HT galaxies- that is, 
the maximum density of the groups containing a HT 
galaxy - show that HT galaxies reside in extremely 
dense regions. Table |3] gives the cut-off densities for 
each of the HT galaxies. We note that due to incom- 
plete optical coverage of our region not all redshifts 
have been determined hence the cut-off densities are 
lower limits. This is especially true for Galaxies A and 
C where there clearly exists close galaxies, but without 
reliable redshifts these galaxies are not incorporated 
in our cut-off density calculation. 

In fact, we find that the HT galaxies are in the 
densest regions of the cluster. We assume that the den- 
sity of the environment has no bearing on whether a 
galaxy is radio-emitting (as per lVenturi et all (|2000l )). 

4.1 Density Cut-OfFs 

To ensure that the extremely dense environments sur- 
rounding our HT galaxies are statistically significant, 
we perform Monte Carlo shuffles on our density cut- 
offs. 
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Table 4. Fraction of Monte Carlo Shuffles that have cut- 
off densities greater than or equal to the mean (or median) 
HT cut-off density. 



Test 




Mean 


Median 


HT V non 


HT radio 


<0.001 


0.001 


HT V non 


HT 


<0.001 


<0.001 


non-radio 


V radio 


0.011 


0.081 



Monte Carlo shuffles work by calculating the 
mean (or median) density cut-off of the five HT galax- 
ies. We then choose five random galaxies and average 
their density cut-offs. This is repeated 100000 times 
and the fraction of shuffles that produces density cut- 
offs greater than the mean (or median) density cut-off 
of the HT galaxies show the statistical significance of 
the HT galaxies' density cut-offs. 

If HT galaxies are truly in denser regions than 
non-HT galaxies we would expect that 

• HT galaxies are found in statistically denser re- 
gions than non-HT galaxies with radio emission, 

• HT galaxies may or may not be found in statis- 
tically denser regions than non-HT galaxies without 
radio emission. This is because galaxies without ra- 
dio emission can potentially be in denser regions than 
HT galaxies but as they have no radio morphology, a 
comparison is meaningless, and 

• radio galaxies are not found in statistically differ - 
ent regions to non-radio galaxies (|Venturi et al.|[2000l ). 

Consequently, we perform Monte Carlo shuffles 
on three different sets of sample populations. 

(i) HT V non-HT radio 

(ii) HT V non-HT 

(iii) radio v non-radio 

The results are shown in Table |4] A result is con- 
sidered statistically significant when the Monte Carlo 
fraction is less than or equal to 0.01. 

The values generated for HT v non-HT radio and 
HT V non-HT are statistically significant for both av- 
erage and median statistics which means that HT 
galaxies do reside in regions of enhanced density. Us- 
ing the mean, it appears that there is a significant dif- 
ference between radio and non-radio galaxies, but this 
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Table 5. Radii of the groups containing HT galaxies. Col- 
umn 1 is the HT galaxy and the remaining columns give 
the projected radius of the group (in kpc) with the specified 
population. 



Table 7. The fraction of Monte Carlo Shuffles that pro- 
duced five groups with an average radius less than the av- 
erage radius of the HT galaxy groups, when all groups are 
considered and when only groups around radio galaxies are 
considered. 



30 



Group Population 



All Radio 



158 
105 
108 



177 

134 
141 



402 
264 
273 
392 
296 



540 
357 
360 
591 



848 
558 
474 
721 



1092 
860 
1172 
1222 
1127 



Table 6. Average radii for all groups. Column 1 is the 
class of galaxy (HT: the five HT galaxies; radio: the 46 
radio galaxies (including the five HT galaxies); all), and 
the remaining columns give the average radii of the groups 
(in kpc) with the specified population. 



2 
3 
4 
5 

10 
15 
20 
30 



0.004 
0.007 
0.005 
0.008 
0.021 
0.073 
0.244 
0.891 



0.045 
0.010 
0.009 
0.010 
0.023 
0.053 
0.160 
0.735 



oup population 
5 10 



49 86 
95 153 
109 163 



211 
216 



154 
261 
268 



325 
504 
489 



473 682 1095 
674 751 1055 
645 749 1020 



significance disappears when using the median statis- 
tic. We attribute this to the very high densities of the 
HT galaxies, which could skew the radio mean, but 
would have little effect on the median. 



From the Monte Carlo Shuffles we can see that 
the smaller radii associated with HT galaxy groups 
is statistically significant for groups with 10 or fewer 
galaxies. There does not appear to be a significant dif- 
ference between the groups around radio galaxies and 
those around all galaxies. The slightly different val- 
ues at the group population of two may be attributed 
to the smaller number of radio galaxies, and thus the 
likelihood of randomly selecting a HT galaxy group is 
higher. 



4.2 Radii 

To verify the results from the density parameter, and 
to explore the environment further out from each of 
the HT galaxies, we conduct an analysis using pro- 
jected radius. For this, the groups are still identified 
using the HIPASS group finder in the usual way, how- 
ever we now examine their projected radii as the num- 
ber of galaxies in the group (group population) in- 
creases. 

Table [5] gives the projected radii of the groups 
found around each HT galaxy, while Table|B]shows the 
average projected radii of the groups found around all 
the galaxies in the catalogue, and is broken up accord- 
ing to galaxy class (ie. HT, radio, all). The average 
group radii for HT galaxy groups is roughly half that 
of those around other classes of galaxies, for groups of 
less than 10 members. This suggests that HT galax- 
ies are surrounded by closer companions than non-HT 
galaxies, and supports the result that HT galaxies re- 
side in denser regions of the cluster. 

We perform Monte Carlo shuffles on our data to 
determine whether the smaller radii associated with 
HT galaxy groups are significant. We calculate the 
average radius of the HT galaxy groups as a func- 
tion of group population, and then choose five ran- 
dom galaxies and calculate their average group radii 
in the same way. This is repeated 100000 times, and 
the fraction of shuffles that produce an average radius 
smaller than the average HT galaxy group radius for 
the same group population is shown in Table [T] 



5 DISCUSSION 

5.1 Close Companions to the HT Galaxies 

Looking more closely at the high density groupings, we 
find that each HT galaxy has a close companion. Using 
DSS images, we locate the closest galaxy as identified 
by the HIPASS group finder; these companion galaxies 
are shown on radio-optical overlays in Figure [S] Each 
HT galaxy has a close companion with projected sep- 
arations ranging from 15 to 119 kpc. 

5.2 Analysis of Statistics 

In our Results section we showed that the densi- 
ties surrounding HT galaxies were significantly higher 
than those surrounding non-HT galaxies. In the fol- 
lowing sections we will discuss the implications of this 
result. 

It has been well documented in the literature that 
there is no correlation between the density of a re- 
gion and the number of radio galaxies observed (e.g. 
IVenturi et al.li2000 ) . Consequently the marginally sig- 
nificant P-values obtained for the non-radio v radio 
mean density test are startling as it suggests that ra- 
dio galaxies exist preferentially in dense environments! 

We are reluctant to believe that radio galaxies 
exist preferentially in dense environments. If anything, 
we expect the densities of the non-HT galaxies to be 
higher than the non-HT radio galaxies as high density 
non-radio galaxies could potentially skew our data. We 
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Figure 5. The HT galaxies and their companion galaxies. The HT galaxy is shown in the "zoomed out" image while the 
optical galaxy and it's companion ar e shown in the "zoomed in" image. The 1.4 GHz ATCA data is red while the 843 
MHz MOST data l|Mauch et aLlbOOal ) is blue. For both frequencies the lowest contour is three times the image RMS noise 
(0.1 mjy beam~^ at 1.4 GHz) and contours increase at \/2 intervals. As noted earlier, the radio galaxy nestled between 
the lobes of Galaxy A is a background source and not considered in this study. The green square identifies the closest 
companion to the HT galaxies as identified by eye. We note that this is the same as the companion identified using the 
HIPASS group finder in three cases. We do not have redshifts for the companion galaxies to Galaxies A and C - this is a 
possible explanation for their lower cut-off densities. The arrow indicates the hypothesized trajectory of the galaxy. 
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attribute this bias to the difFerence in optical and radio 
coverage. The optical data have an average coverage 
of ~90% down to a magnitude of 18. Conversely, our 
radio coverage is complete to flux densities of > 0.5 
mjy at 1.4 GHz. This is essentially 100% complete 
for HT galaxies as the faintest HT galaxy has a flux 
density of 11 mJy at 1.4 GHz. The incomplete optical 
catalogue has the consequence of skewing the density 
of the optical catalogue downwards; that is, it appears 
less dense. 

From Figure [1] we can see that the optical cov- 
erage is 100% down to mi,=18 over most of the field 
and the average coverage decreases uniformly over the 
field with fainter magnitude cutoffs. This means that 
the optical coverage surrounding the HT galaxies is 
not any different to the optical coverage surrounding 
the other radio galaxies. Thus, we think that the HT 
V non-HT test is free of optical biases. 

The p-values obtained for both HT galaxy tests 
strongly suggest that the densities surrounding HT 
galaxies are significant and hence we are confident (to 
~99%) that HT galaxies exist in significantly denser 
environments than non-HT galaxies. 

The results derived from the 2D density enhance- 
ment analysis using projected radii appear to support 
the results we obtain from our 3D density enhance- 
ment study. 

Our results also show that the tighter-knit envi- 
ronments surrounding HT galaxies remain significant 
for group populations of up to 10 galaxies. This sug- 
gests that high densities are required for HT galaxy 
morphology to occur. 

5.3 A Brief Foray into the Shapley 
Supercluster 

To ensure the density values obtained for the HT 
galaxies in the central region of the HRS are realis- 
tic, we test our process on the closest supercluster in 
the local universe: the Shapley supercluster. We do 
not have data for th is region so w e rel y on the radio 
data available from iMillerl (|2005l ) and IVenturi et al.l 
(|2000l ). From their papers we identify five HT galax- 
ies in the three central clusters of the Shapley super- 
cluster. (The three clusters are a well-studied merger 
A3556-A3558-A3562, not dissimilar to A3125/A3128 
in the HRS.) The positional and redshift catalogue 
is taken straight from NED, so our visual catalogue 
compilation method is not utilized. The optical cover- 
age of the Shapley supercluster is com parable to our 
data for the HRS jVenturi et al.ll2000l '). Furthermore, 
the sensitivity of the radio data is also comparable to 
our radio data (jVenturi et al.| [2"o00l) so we will assume 
the two superclusters may be directly compared. The 
Shapley supercluste r has roughly dou ble the galaxy 
density of the HRS (|Proust et aLlbOOel ). 

The density cut-off values for the HT galaxies in 
the Shapley supercluster are similar to the values ob- 
tained for the HT galaxies in the HRS (Table [HI . 

The values for the Shapley supercluster are 
slightly higher than for the HRS but this may be at- 
tributed to the denser environment of the Shapley su- 



Table 8. Cut-ofT densities for HT galaxies in Shapley. 



Vi = 1000 km s- 
HT galaxy Cut-ofT density 

(galaxies /Mpc^ ) 

J132357-313845 8639 

J133331-314058 1708 

J132802-314521 183654 

J133542-315354 367 

J133048-314325 16729 



percluster. However, the similarity of the results for a 
peculiar velocity of 1000 km s^^ in both the HRS and 
the Shapley supercluster gives us confidence in the ve- 
racity of our results. Thus we may say with confidence 
that HT galaxies do form in regions of enhanced den- 
sity. 

5.4 Further Evidence for High-Density 
Environments 

Inspection of archival ROSAT X-ray data shows that 
four of the five HT galaxies (A to D) reside in regions 
of enhanced emission near the centres of A3125 and 
A3128. There is also evidence that the fifth galaxy 
(galaxy E) is in a region of X-ray emission which is 
above the X-ray background but it is at a much lower 
level than for the other galaxies. Since X-ray emission 
also traces regions of enhanced density these results 
support th e statistical analysi s presented above. In- 
terestingly, iBurns et al.l (|l994l ) found a high correla- 
tion between X-ray emission and radio galaxies, and 
that X-ray emission aro und WATs was g e nerall y more 
extended. Furthermore, I Miller fc OwenI l|l999l ) found 
that X-ray emission around FRI galaxies is more ex- 
tended than that around FRIIs. Our results are consis- 
tent with these findings. Details of the X-ray emission 
will be discussed in a separate paper. 

The morphology of the HT galaxies is itself in- 
dicative of a dense environment, as it is believed that 
the jets of FRI galaxies are less efficient in transport- 
ing energy due to the dense ICM surrounding them 
l|Fanaroff fc Rilevlll974l ). 

The morphology of HT galaxies has always been 
believed to be due to either: 

• the host galax y having a high p eculiar velocity 
(radio trail theory, (iMilev et al.lll972l )'). or 

• the host gal axy being s wept by strong winds 
(cluster weather, (|Burnsill998i ')'). 

Using the optical data for our HT galaxies we 
find that the velocities are not exceptionally different 
to the cluster mean and hence it would appear that 
the host galaxies are victims of cluster weather. How- 
ever, studies of cluster weather indicate that if cluster 
winds were the sole cause of the HT morphology, the 
velocity requir ed for ram pre ssure is less than previ- 
ously thought l|Pinknevlll995l ). While we find that the 
velocities of the HT galaxies are not especially great. 
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the densities in which they reside are among the dens- 
est regions of the cluster. Although we measure the 
number density of galaxies, this value follows from the 
cluster potential, and thus should also be representa- 
tive of the density of the ICM due to gravitational 
attraction. 

We propose that the morphology is due to ram 
pressure effects, but it is the exceptionally high gas 
densities that contribute to the morphology, as op- 
posed to peculiar velocities or high velocity winds due 
to cluster weather (as Pramocpf^). 



6 CONCLUSIONS 

We have compiled the largest positional and spectral 
catalogue of the major merger A3125/A3128. This 
well studied merger lies in the central region of the 
HRS, the second largest supercluster in the local uni- 
verse. Using 1.4 GHz radio data we identified five HT 
galaxies in this region. 

To investigate the relationship between the HT 
galaxies and their underlying dynamical structure we 
identified substructure in the central region of the 
HRS. We found, using the HIPASS group finder, that 
the HT galaxies all exist in some of the densest re- 
gions of the cluster merger system A3125/A3128. We 
performed a series of statistical tests to probe three 
populations of galaxies within A3125/A3128 (15,000 
kms"^ < cz < 35,000 kms"^): all galaxies (582), all 
radio galaxies with a spectroscopic redshift within 
A3125/A3128 (46) and all HT galaxies (5). Results 
of the single-tailed permutation tests showed, to 99% 
confidence, that HT galaxies reside in significantly 
denser environments than either non-HT galaxies or 
radio galaxies which are known members of the clus- 
ters. These results are also replicated in a comparative 
study of the five HT galaxies in the centre of the Shap- 
ley Supercluster. 

Additionally, we found each of the HT galaxies 
in A3125/A3128 reside in regions of enhanced X-ray 
emission and have a very close companion (ranging 
from 15 - 119 kpc in projected separation). As the 
peculiar velocities of the host galaxies are not signifi- 
cantly different to the mean cluster velocity this sug- 
gests that it is the density of the environment which 
plays a significant part in the bending of the jets. How- 
ever, the presence of the close companions is also likely 
to play a role in at least the motion of the host galax- 
ies, though it is likely that the morphology of the jets 
on kpc-scales is more closely related to ram pressure 
and hence cluster weather. Thus, we propose that the 
combination of both a close companion and the clus- 
ter environment contribute to the bent morphology of 
HT galaxies. 

Ultimately, we have found that HT galaxies are 
excellent indicators of dense regions and dynamical 
interactions within galaxy clusters and thus are clear 
signposts of "cluster weather". As a result, HT galax- 
ies may be used as barometers for the storms that 
occur in the cluster environment, although due to 
the complication of close companion galaxies, reading 



these barometers may not be as easy as first thought. 
We propose that the diffuse tails, on at least kpc- 
scales, can be used as giant wind socks and are ef- 
fective beacons of cluster weather. Detection of the 
full extent of HT galaxies requires low frequency ob- 
servations of reasonable resolution. Instruments such 
as the GMRT and LOFAR will provide large samples 
of HT galaxies which may then be used as beacons of 
high-density regions in clusters. 
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